Theory of Non-Coherent Spin Pumps 
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We study electron pumps in the absence of interference effects paying attention to the spin de- 
gree of freedom. Electron-electron exchange interactions combined with a variation of external 
parameters, such as magnetic field and gate potentials, affect the compressibility-spin-tensor whose 
components determine the non-coherent parts of the charge and spin pumped currents. An appro- 
priate choice of the trajectory in the parameter space generates an arbitrary ratio of spin to charge 
pumped currents. After showing that the addition of dephasing leads to a full quantum coherent 
system diminishes the interference contribution, but leaves the non-coherent (classical) contribu- 
tion intact, we apply the theory of the classical term for several examples. We show that when 
exchange interactions are included one can construct a source of pure spin current, with a constant 
magnetic field and a periodic variation of gate potentials only. We discuss the possibility to observe 
it experimentally in GaAs heterostructures. 

PACS numbers: 72.25.-b, 73.23.-b 



I. INTRODUCTION 

A pump is a very common device, it appears in many 
shapes and types, in physical systems as well as biolog- 
ical systems. A pump converts a periodic variation of 
the parameters that control it into a direct current. For 
example, in an Archimedes (or Auger) screw pump after 
one revolution "buckets" of water are lifted up, while the 
blades of the screw return to their position at the initial 
stage of the revolution^ 

As modern electronic devices, including electron 
pumps, become smaller and get into the mesoscopic size 
of a few nanometer to a fraction of micrometer, we need 
to consider also the wavy, i.e., quantum mechanical, na- 
ture of the electrons. Similar to the classical devices, a 
quantum pump converts periodic variations of its param- 
eters at frequency u> into a DC current^ 

In both quantum and classical devices the pumping 
rate is proportional to the liquid pumped in one turn 
of the pump and to the revolution rate lo. When the 
Archimedean screw is rotated too fast, turbulence and 
sloshing prevents the buckets from being filled and the 
pump stops to operate. Similarly, in small electronic de- 
vices, the pumping rate is proportional to u> only in the 
adiabatic limit — when u> is small enough. 

While pumps exist for both interacting coherent quan- 
tum systems and interacting classical systems, the main 
theoretical studies of mesoscopic pumps were concen- 
trated on the wavy nature of the electrons describ- 
ing the pumps in terms of non-interacting coherent- 
scattering theoryiiii* Several works discuss the effect 
of dephasingiS** which smears the wavy nature of elec- 
trons, suppresses interference effects and renders the sys- 
tem non-coherent, i.e., classical. The complexity of the 
full quantum problem in presence of interactions*^ al- 
lowed its study only in few examples of open quantum 
dotsii and Luttinger liquidii* 

On the other hand, the early experimental studies of 
quantum-dot-turnstile pumpA*** are described in clas- 
sical terms of interacting systems, namely, capacitance 



and oscillating "resistances" of tunnel barriers. In a 
more recent experiment ^* coherent pumping was ob- 
served. (However, parasitical rectification effects may be 
relevant i**i*) 

In this manuscript we study the classical non-coherent 
effect of pumping and in particular we show how this clas- 
sical effect directly emerges out of the quantum mechan- 
ical formulation* when dephasing sources exist. While 
the quantum-scattering description may include physics 
related to the spin degree of freedomji* it treats electron- 
electron interaction on the Hartree level only. In this 
manuscript we include both direct and exchange inter- 
action within the framework of a classical theory. This 
enables us to suggest a scheme to build a pure "spin bat- 
tery" (see Sec. IV C|l which is a key concept in the field of 
spintronics. 

The remainder of the paper is organized as follows: 
in Sec. [H] we derive an expression for the non-coherent 
contribution for small electron pumps using a classical 
model that neglects interference effects completely. The 
formulation of pumps in classical terms allows us to in- 
clude rather easily the effects of interactions between the 
electrons. 

To describe effects of dephasing in a controlled man- 
ner**» we include in Sec. IIIII following Ref. 0, the effect 
of additional voltage leads on the (non-interacting adi- 
abatic) quantum scattering theory of pumpsi*i*»i* We 
show that when the coupling to the dephasing leads is 
tuned properly to cause complete dephasing, the Brouwer 
formula** which relates the pumping current to the scat- 
tering matrix and its time derivatives, reduces to the clas- 
sical expressions developed in Sec. |n] We compare the 
magnitudes of the quantum and classical contributions 
and study under what conditions the classical contribu- 
tion dominates. 

After constructing and justifying the classical theory of 
non-coherent pumps, we generalize it in Sec. IIVI to deal 
with spinfull electrons endowing our result with a topo- 
logical interpretation. Finally, in Sec. we deal with 
a possible realization of spin pumps in two dimensional 
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electron gas of GaAs heterostructures. We show that 
under rather general conditions the spin current (simi- 
lar to the Einstein's DC conductance formula) is given 
in terms of the thermodynamic density of states tensor 
of the system (the compressibility tensor). In the pres- 
ence of a constant magnetic field, we find an appropriate 
trajectory in the parameter space such that a pure spin 
current will flow. This effect vanishes in the absence of 
exchange interaction. 

In appendixlAlwe discuss the relation between the clas- 
sical non-coherent pumped current and the biased cur- 
rent, generated by rectification effects^ In appendix [B] 
we explore the relation between our theory and the the- 
ory of non-linear response. In this manuscript we do not 
include effects of charge discreteness leaving this for a 
future study>i^ 



II. CLASSICAL DESCRIPTION OF 
NON-COHERENT PUMPS 

Consider the electrical circuit depicted in Fig. ^ Its 
purpose is to charge the capacitor C from left and then 
to discharge it to the right, thus producing a net current 
from left to right without a bias voltage. The gate volt- 
age, V g , periodically charges and discharges the capacitor 
while the resistors, Rl and Rr, control the direction of 
the charging and discharging processes. 

To analyze the system it is convenient to define an 
asymmetry parameter 



a{R L ,R R ) = \ 



Ri 



2 Rr + Rl 



(1) 



running between 1/2 for Rr ^> Rl and —1/2 for Rr <C 
Rl- 



R L(t) | L 




FIG. 1: Equivalent circuit of a classical pump. A typical 
pumping cycle occurs when we charge the capacitor through 
the left resistor (keeping Rr ^> Rl) and discharge it through 
the right resistor (when Rr <C -Rl). In that case a net DC 
current flows in the clockwise direction in the large current 
loop. 



The pumping circuit is operational even at a con- 
stant gate voltage, when V g = Vq. Suppose that the 



capacitor may be typified by a parallel plate capacitor 
with a tunable area A, separation d and permittivity 
e. Initially the capacitor is at equilibrium with charge 
Qo = CV = ^f. Then its area is varied to A' = A- 5 A 
while a « — 1/2 (Rr -C Rl). As a result a charge 5Qo 
will leave the capacitor to the right until equilibrium is 
restored. Then by changing the area back from A' to A 
when a f=a 1/2, we will pull the same amount of charge 
SQo into the capacitor from the left. Repeating this pro- 
cess with a period r will result in an average pumped 
current SCVq/t flowing from left to rights (A periodic 
variation of d or e will have a similar effect.) 

Proceeding formally, let Q(t) be the instantaneous 
charge on the capacitor. The charge leaving the capaci- 
tor during a small time interval dt is —dQ(t) = —Q(t)dt. 
From current conservation and Kirchhoff's rules: 



I L + Ir = Q(t) = Ic 



IlRl — IrRr, 



(2) 



the fraction of this out going charge, —dQ(t), leaving via 
Rl or via Rr is 



dQ R {L){t) = - 



R 



L(R) 



Rl + Rr 



dQ(t). 



_ (Il)-(Ir) 



(3) 



where 



We define the pumped current as I 
(O) = r _1 JZ 0(i)dt denotes average over a period. Then 
/ may be expressed in terms of Q(t) as 



I= l - 

T 



dt-(h 



LR = - 



dta(t)Q(t). (4) 



The charge Q(t) should be determined by the equa- 
tions of motion of the system which are governed by 
a Lagrangian C, including a source term in the Euler- 
Lagrange equations which introduces dissipation: 



SC(Q,Q) dS£(Q,Q) 



SQ 



dt 



= RQ, 



(•5) 



where R = R L \\ Rr = R l Rr/(Rl + Rr)- 

If the parameters, x(t),y(t), . . . , that control Q and 
a in Eq. Q are varied slow enough then Q and a are 
functions of the instantaneous value of these parameters, 
and do not depend explicitly on time: 



Q(t) = Q[x(t),y(t), 
a(t) = a[x(t),y(t), 



(6) 



The parameters x(t), y(t), . . . can be for example 
V g , d, A, e, Rl, Rr or any combination of them, e.g., a 
itself. We will refer to this slow limit as the adiabatic 
limit. For each case that we study we will check how 
large should be the period r for the adiabatic limit to be 
established. Roughly, the adiabatic condition is estab- 
lished when t is larger than the effective RC time of the 
circuit. 

If there are only two parameters, then in the adiabatic 
limit the pumped current is 



1 

I = - 



dta[x(t),y(t)]Q[x(t),y{t)}. 



(7) 



3 




FIG. 2: The pumped charge per period is the flux of an effec- 
tive magnetic field inside the parameters trajectory. 



Using now 



dQ dQ dx dQ dy 
dt dx dt dy dt ' 



the current can be rewritten as a line integral along a 
trajectory L in the parameter space (see Fig. EJ: 



(ad x Q, adyQ) ■ (dx, dy). 



(8) 



Because the parameters are varied periodically in time 
the trajectory is closed. Using Stoke's theorem, the line 
integral can be transformed into a surface integral on the 
closed surface S bounded by the trajectory L, 



I = 
V x A cB = 
A cS = 



1 
r 



dxdyB cS = 



xcff 



d x ad y Q — dyCid x Q = Va x VQ, 



aVQ or - (aVQ 



QVa), 



(9) 



where B cS is an effective magnetic field in parameter 
space. The last ambiguity in the definition of A eS is a 
result of a gauge freedom: an addition of the gradient 
term iV(aQ) does not change B cS . 

Equation © suggests that the charge pumped per cy- 
cle, rJ, is equal to the flux of the effective magnetic field 
through the loop in parameter space <f> , as depicted in 
Fig. [5] Similar topological formulation for the pumped 
current was discussed in the quantum case^ while here 
we obtained a similar structure for the classical situation. 



A. Example: x = a, y = V g 

Since no current is pumped if the asymmetry param- 
eter a is kept fixed, we consider the simplest example 
where a itself is a pumping parameter. Substituting the 



Lagrangian of the circuit depicted in Fig. ^ (generalized 
to the case where the capacitance depends on its charge 
or on the voltage across it): 



£(Q,Q) = -E C (Q) + QV, 



rQ 

E C (Q) = / dQ'V c (Q') 
Jo 

dQ' 



VrXQ) 



o C(Q') 



(10) 



in Eq. one gets the equation of motion V g — V C (Q) 
QR. Using C — dQ/dV c we get 



Va-V P = RCVr,. 



(11) 



The meaning of the adiabatic limit is clear in this equa- 
tion. The typical time scale for changes in V c is the pe- 
riod time r. Therefore if RC <C t we can neglect the 
right hand side of the equation and establish the adia- 
batic limit V c (t) = V g (t). Substituting this into Eq. © 
one gets 



i = \ H rn>w/u ; 



(12) 



so that C plays the role of the effective magnetic field. 

For C = 10~ 15 F at a frequency t _1 = 1GHz and with 
a gate voltage oscillations of lmV the maximal pumped 
current is In A. 



B. Application to a more general electrical circuit 

Consider the electrical circuit shown in Fig.[3]that gen- 
eralizes the circuit depicted in Fig. ^ The analysis of the 
two circuits is similar: let us assume that we change the 
charge on capacitor C mo by dQ mo (by changing V g mo ), 
while keeping the charge on the other capacitors con- 
stant. Using Kirchhoff's rules one can show that the 
fraction of dQ mo flowing to the left via R\, dQ^, or to 
the right via Rn, dQn, is: 



dQ R {L){t) 



r>L(R) 
Em=l R > 



-dQ mo (t), 



(13) 



where R^ o = YZ°=i Rm and R* = £m=m +i R m are 
the resistances to the left and to the right of capacitor 
C mo respectively. 

The superposition principle generalizes Eq. to 



dQ R { L ){t) 



E k R L k iR) dQk(t) 

Z^m=l Jx m 



(14) 



for any variation of {V g k}- Similarly, Eq. (|12fl which 
holds in the adiabatic limit becomes 



N-i 

- // V C(V g k )da k dV g fc , 
T JJs k=i 



(15) 
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where ctk = hjrkz 



In the next section we will show that a small quantum 
system subjected to dephasing can be described by the 
circuit depicted in Fig. |3J 



> i > 



v 01 



: Cn.i 

- V.N-, 



III. 



FIG. 3: Generalization of the circuit in Fig. 



CONNECTION WITH QUANTUM PUMPS 
SUBJECTED TO DEPHASING 



The scattering approach for pumps 4 - relates the DC 
current flowing through a time dependent scatterer with 
its scattering matrix and with the time derivatives of 
the scattering matrix. However, in real physical sys- 
tems there are processes that lead to uncertainty in the 
phase of the single electron; for example interaction with 
phonons and interaction between the electrons mediated 
by the electromagnetic environment. This effect is called 
dephasing. It is expected that a classical description will 
gradually take place as the dephasing becomes stronger. 

To study the effect of dephasing on the quantum coher- 
ent pump in a controllable way we introduce dephasing 
leads&SSi and show that in the limit of strong dephasing 
the pumping current approaches the classical limit given 
in Eq. Q15[l. The effect of dephasing can be described 
in different ways which differ in details, but this variety 
of ways does not alter the main conclusion that a quan- 
tum system with strong dephasing can be described by a 
classical theory. 

Consider a conducting wire subjected to a gate poten- 
tial V g (r), r e [0, L], as shown in Fig. 0Ja). To introduce 
dephasing in a controlled way we connect N — 1 wave 
splitters of 4 legs2i at the points r = i£, i = 1,2, ...,N — 1 
along the wire, as shown in Fig. Efb) (for N = 3). The 
length £ determines the dephasing length of the model 
together with the wave splitter parameter e, as we will 
explain later. The wave splitters are described by the 
scattering matrix 



Ssplitter(e) 



/ VT 



V o 










o 



where the third and forth lines and columns of the matrix 
correspond to the two legs of the reservoir that serves 
as a dephasor. Each wave splitter i is connected to a 
reservoir i held at voltage V.. To mimic a dephasor that 
influences only the phase coherence of the waves in the 
sample and does not influence the total current, we tune 



the voltages {Vi} so that no net charge flows into any of 
the reservoirs. In our time dependent problem we will 
assume that such conditions hold at all times. 

(In a similar model of dephasing 2 ^ one introduces an 
extra phase </> m selected points, which is averaged out 
after squaring the desired amplitude. The approach that 
we use differs from the later in the following way: the 
introduction of the reservoirs is accompanied by the ad- 
dition of contact resistances, which renders the total re- 
sistance from left to right higher.) 



V q (r) 



(a) 



wave 
v splitter 



left Va=0 


L 




\S(E) 








R 


right VyO 
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(b) 
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FIG. 4: (a) The gate potential along the original quantum 
system without dephasing. (b) Introduction of dephasors, 
controlled by a parameter e. (c) e = 0: The system is totaly 
coherent, (d) e = 1: This configuration doesn't allow any 
interference effects between different scatterers S m - 



For general e the whole system, including the dephas- 
ing leads, is described by a S matrix of rank 2N. The 
current in each lead is composed of two contributions, the 
pumped current (Brouwer formula 4 ") and the biased cur- 
rent (Landauer-Buttiker-Imry conductance formula 23 ). 
The voltages {Vi} produce biased currents which cancel 
the pumped currents in each reservoir. 

In the limit e — > (without dephaing) the S matrix is 
block diagonal: it is a direct sum of N rank 2 matrices. 
One of them connects between the left and right reser- 
voirs and describes scattering from the potential V g (r) 
between r = and r — L. The other matrices trivially 
connect each reservoir only to itself, i.e., the reservoirs 
are not connected to the wire, see Fig. IHc). 

In the limit e — > 1 we again have S = ©„ =1 S rra . Now 
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S m describes a potential barrier V g < m ) ( r ) between two 
dephasors given by 



aim) 



M -fV g (r) (m-l)Kr< 
[ ' i else 



ml 



(16) 



Solving Eq. (29 for we find>2i 

, n _ (1/2 + a-m a )Rm + Engmp+1 j^ro j^. 



The original barrier is split into N sections, see 
Fig. Hid). Notice that while the phase coherence be- 
tween different sections is lost, they are correlated by 
the requirement for zero current in the reservoirs. The 
dephasing length is thus equal to I. 

Let us concentrate on the e = 1 case which is easier to 
handle, (mentioning that, in principle the model could be 
solved for any e) and parameterize the unitary scattering 
matrix of each section as 



2wm=l IV - 

aQR = — ff — —q aQm - (22) 

Lm=l 

We thus obtain Eq. © with R L -> YZ^i ^+(V 2 ~ 



* mo )R mQ and R R -> (1/2 + aQjRQ c 



m— mo + 



R Q 

1 -"-m* 



,»Am [ el??m cos ®m ie l< ^ m sin # m 
I i e t <t>™ sin6* m e~ lrim cos# m 



For this parametrization the Brouwer formula reads^ 

27r/e dQ m ^ L = -dA m - cos 2 9 m d-q m + sin 2 6 m d<j) m 
2n/e dQ mt R = -dk m + cos 2 9 rn di] m - sin 2 m d<f> m , (18) 

where dQ m ,L and dQ m .R are the pumped charges flowing 
out of section to to the left and to the right respectively, 
see Fig. 0] 

To understand how pumping takes place in the e = 1 
case suppose first we change V g (r) only in section too 
such that only the parameters A mo , 9 mo , f]m and 4>m a 
may change. Consequently there are pumped charges 
flowing out of (or into) section Too to the left and to the 
right in leads mo,L and too, R-, given by Eq. (|18fl . The 
change in the charge of section too (by analogy with the 
charge of the capacitor in Sec. ITT|> . is 



During a generic pumping process the potential is var- 
ied in different places, i.e., we have to consider simulta- 
neous variations of V g (r) in sections m ^ too. Since the 
different sections are connected classically when e = 1, 
(17) and Eq. (121 H is linear in the source term oc Qm , we sum 
over all the contributions to dQh and dQx arising from 
each section. We then obtain Eq. (|14H in which the re- 
sistors are obtained from the Landauer resistances {R m } 
and the quantum partitioning coefficients {a^} as 



R m = (1/2 + < 



)R W 



(1/2- a«)i& (23) 



dQ r . 



-\dQ mo ,L + dQ mo> n). 



(19) 



We denote its partitioning into left and right by —(1/2- 



dQ r , 



dQ mo ,L and -(1/2 - a;« ) dQ mo = that for this case dQ L 



First, notice that for large N Eq. (1221) is insensitive to 
the actual value of {a%}. It occurs because a change in 
these values modifies only one term in the sum. Second, 
as we already mentioned above, this mapping is impossi- 
ble when dQ m „.L — —dQ m0t R, i.e., when the effect of the 
variations in the potential in a single coherent section is 
to transfer charge from left to right, without charging nor 
discharging this section. This situation occurs only for a 
specially tuned asymmetric potential inside a single co- 
herent section. However, our equivalent circuit requires 
a nearly constant potential in each section, as there is 
only one gate in each section (see Fig.©. Notice that for 
large N this is a small effect: one can see from Eq. (|22f) 

R$ 



dQm ,Ri where a m is called the quantum partitioning 



El 



-dQ 



rriQ,L(R) • 



coefficient of section m . (When dQ r . 



one should 



write the formulae below explicitly in terms of dQ mo .L 
and dQ mo ,R-) 

The voltages on the reservoirs to = 1, . . . , N — 1 are 
adjusted to cancel out the pumping contribution, so that 
there is zero current in each reservoir. That gives the 
equations 

V m — V m -i V m — V m +\ 



R 



Q 



R Q 



-(5 m , mo _i(l/2 + a« ) + 5 m , mo (l/2-a£ ))Q mo ; 

to = 1,2,..., TV - 1, (20) 

where R m = (tt sin 2 ^m) 1 are the quantum resistances 
of the different sections. The currents into the left and 
the right reservoirs are 



Vx/R?dt - S moil (1/2 + ag)dQ mo , 



A. Comparison between coherent and 
non-coherent effects 



The dephasors connect incoherently (or classically) 
sections of length £, which conserve internal phase coher- 
ence. Indeed Eq. (1221 is derived by classical circuit the- 
ory, however, coherence effects still determine its param- 
eters (dQ mo and the resistances {R m })- In this section 
we separate coherent from non-coherent contributions to 
dQm and compare between their magnitudes. 

To do so, suppose the gates length i gato 3> I such 
that the characteristic length scale for changes of V g (r) 
is larger than i. Then V g (r) can be approximated by 
Eq. I|16|) where V g i m ) (f ) is a rectangular barrier of width 
I and height V g[m) = V g (£{m - 1/2)). The transmission 
coefficient of section m is thus 



dQ R = V N -i/R%dt - b ma . N (1/2 -a mo ) dQr. 



(21) 



ik,„ e 



(fci 



(24) 
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where kp is the Fermi wave number and k m = 
y/2m(Ep — V g ( m ))/h 2 is the wave number inside the bar- 
rier of section m. (Since all the potential voltages oc uj, 
for very small u> the Fermi wave numbers in all sections 
are similar and we ignore the difference between them.) 

Consider one of the sections that are influenced by the 
long gate, say mo- Due to the inversion symmetry of each 



barrier we have r/ n 



and 



0, and thus A, 



arg(i mo ). As V g i mo ) increases, particles are pushed out 
symmetrically (in this case a^ = 0), and by Eq. Ijl8(l 



their charge is dQ mo ,L+dQ mo ,R 



dQ 



m 



-edA mo /n 



Eq. (|24[1 can be derived by summing over all possible wave 
trajectories from the left to the right of the barrier; the 
numerator corresponds to the straight (classical) path 
through the barrier while the denominator appears after 
summing the rest of the paths. Following this observation 
we artificially decompose dQ ma into a non-coherent part 
[numerator of Eq. 124fl ] and coherent part [denominator 
of Eq. (EU] : 



dQ, 



edV g ( mo ) 
where J) m , 



-6{E f 



d arg /„ 



irdV, 



g(m ) 



(25) 



is the density of states at wave 
number k mo and f mo is the denominator of t rno in 
Eq. (|24|) . We want to compare these two terms for a 
finite change in gate voltage V g ( mo y 

Since \kp + k m \ > \kp — k m \ it is straight forward to 
realize that arg/ mo is an oscillatory function of V^( mo ) 
whose amplitude of oscillations is smaller than ir. There- 
fore the contribution of the second term in Eq. (|25ll . a 
consequence of interference of many paths, is bounded by 
the single electron charge, e, for any variation of V g ( mo y 
Now we have to sum over the quantum contributions 
of the ig a to/^ sections that are influenced by the gate. 
Due to the oscillatory nature of this term, if the gate 
length Zg a t c ^S> £ and we assume the oscillations to be 
uncorrelated then the contribution of the coherent term 
is 0(\/U^e). 

On the other hand the first classical term, which is 
simply the number of states below the Fermi level in a 
ID box of size £ with potential V^( mo ), is not bounded. 
Furthermore the contributions due to various sections of 
length £ which are subjected to the influence of the same 
gate potential add up and give charge 0(l SB ±e/£ qe) where 
q is the average number of states added to a section of 
length £ due to the change in the gate potential. 

We see that the non-coherent (classical) effect domi- 
nates when Zg a to S> £, or when q 1. In that situa- 
tion coherent effects are unimportant and our electrical 
treatment of pumping is relevant; the pumped current 
is then given by minus the expression in Eq. I|15|) with 



C{V) = 



et^Jrri 



-9(E F - eV). 



Trhy/2(E F -eV) 

On the other hand, the classical term can be tuned 
to zero if we oscillate several gates together, such that 
no net charge flows into or out of the system, and then 
the current results only from quantum interference. This 
happens for example when two nearby gates change in 



opposite directions, or in the two dimensional case when 
the gates change the shape of the capacitor but not its 
area. 

To conclude, a quantum pump with dephasors can be 
described by the electrical circuit of Fig.|3where each sec- 
tion of the circuit describes a coherent section of length 
£ of the quantum system. The components of the circuit 
(resistors, capacitors, etc.) are generally determined by 
coherent effects (for example a capacitor Ci might de- 
pend on the voltage V g i in an oscillatory way). In con- 
trast to coherent systems, for the classical circuit (i) the 
charge inside a coherent section is determined only by the 
potential in that region, Eqs. I|18[l and l|19l) ; (ii) the left- 
right partitioning of an extra charge flowing out of (or 
into) a coherent section is simply the classical partition 
of current through two parallel resistors, Eq. i|22[l : (iii) 
the superposition principle holds for arbitrary change of 
the parameters in different sections. 



B. Example 

To illustrate the difference between the classical and 
quantum aspects of pumping consider as an example the 
potential V(r) = 7<5(r) + U9[r(L — r)]. For a square 
shaped path of the parameters (7, U) passing through the 
points (0,0), (oo,0), (00, 5U) and (0, <5f7), the pumped 
current is4 



/ = 



eLuSU 
8n 2 k F 
0{5U 2 ) 



jSU 



\%ir 2 k 2 F 



[nsm 2 (k F L) - sin(2k F L)) 



(26) 



t2 

with ^r = l. Let us add N — 1 dephasing leads as de- 
scribed above and tune e = 1. The matrix Si describes 
the original barrier shrinked to length £ instead of L while 
the matrices S m , m = 2,...,N describe rectangular bar- 
riers of height U and width £. To find the current in the 
presence of dephasing we have to follow the 4 stages of 
the period using Eq. I|22|) and Eq. (|18|l . For example, 
in the third part, where the parameters (7, U) go from 
(00, 5U) to (0, 6U), only the phases Ai, 81, and 771 change. 
To first order in SU, the resistances of the N — 1 barriers 
are equal to the quantum resistance, and thus 



sec 2 9i 



sec 2 Oi+N-1 
(^1 - cos 2 01— ) — 

C?7 G?7 47T ' 



(27) 



where 9, 7/, A and 77 are defined according to Eq. \17}. 
Performing similar calculations for the rest of the period, 
one obtains a DC current, 



eLuSU etudU ,n sin 2 (kp£) 

87T 2 kp 

(D(6U 2 ) 



16ir 2 k 2 F 



N 



- sin(2M0(2 - 1)) 
(28) 



directed to the right. 
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Let us compare Eq. l(2l)|) with Eq. l(2*5|) . The first terms 
are identical, independently of N. These terms are the 
classical contribution which is unaffected by the dephas- 
ing, being related only to local density of states. The 
second terms coincide for N = 1 but differ otherwise. 
The factor sin(2fcL) in Eq. (|23> was transformed into 
sin(2fc£), i.e., coherent effects are restricted to the de- 
phasing length. [This result depends on the dephasors 
position: a single dephasor at r — + is enough to de- 
stroy the interference term completely, since to 0(5U) 
the only section giving rise to interference effects is the 
one with the <5-function.] 

The ratio between the classical and the interference 
terms is of the order of L/Xf, confirming the importance 
of the classical term when L ^S> \f (or Z ga tc 3> Xf )• 

There are both interference terms 0(N ) which survive 
the limit N -> oo, and others, ©(A^ 1 ) and 0(N' 1/2 ), 
which disappear in that limit. The reason for having 
interference effects in the large N limit in this example is 
the presence of the infinite barrier: when the parameters 
(7, U) go from (00, 0) to (00, SU), all the charge repelled 
from sections m = 1, ...,N — 1 is driven to the right. 
Among all these sections, only section m = 1 gives rise 
to the interference term which is independent of N. On 
the other hand when the parameters move from (00, SU) 
to (0, SU), SQr is given by Eq. (|27[) and we see that for 
iV = oowc get Iji = 0. 

IV. SPIN POLARIZED DC CURRENT 

In Sec. [H] we analyzed the time evolution of an elec- 
trical circuit and derived an expression for the pumping 
of charge in the adiabatic limit. In Sec. If f f I we used a 
dephasing model to show that in the limit of strong de- 
phasing the non interacting coherent pumping expression 
is reduced to the classical expression. In this section we 
will generalize the classical equations to include spin. We 
will assume (i) the dephasing length is smaller than the 
size of the system, I <C L, so that the classical expressions 
are valid; (ii) the spin is conserved during a pumping cy- 
cle, i.e., t <§C T s f , with r s f being the mean spin-flip time. 

The generalization of the spinless treatment of Sec. [D] 
to the spinfull case is done by introducing a Lagrangian 
that depends on charges with spin up and spin down. 
The equations of motion with the dissipation term are 
[cf. Eq. ©] 

SLjQaM d 8C{Q a ,Q a ) A 

In 7j7 7a = ttQ<j (29) 

SQcr at SQ a 

where we have assumed that R doesn't depend on spin. 
Now the current and voltages have two components, re- 
ferring to spin up and down. The analysis goes in paral- 
lel to the spinless case: given the time dependent charge 
with spin index a on the capacitor, Q a (t), the pumped 
current is 

I a = - [ dta(t)Q a {t). (30) 



In the adiabatic limit the spinfull form of Eq. is 
^ = lJJ s d ^-Bf, Bf = VaxVQ rT , 

/c ' S = rjj s d ^' ^ ' ^ = ^ ± ^ ' (31) 

where 7 C , S = ij±7|. The 2 dimensional integrals are done 
on the area § bounded by the trajectory L in the param- 
eter space. The V symbol denotes partial derivative with 
respect to the parameters (2, y, z), V = (d x ,d y , d z ). The 
increase of the number of parameters from 2 to 3, cf. 
Eq. J3J, reflects the need to control the spin polarization 
of the pumped current. 

To understand better the meaning of the effective mag- 
netic field let us decompose it into field lines, as 
shown schematically in Fig. [5] According to Eq. I|3I() 
the total charge of spin a pumped from left to right per 
cycle, associated with a given loop L, is the flux of field 
lines of spin a inside the loop L. To illustrate this point, 
the loops Lj and Lj in Fig. EJb) correspond to pumping 
of only spin up and spin down, respectively. 

In Figs. |3Jc) and OJd) we show the charge and spin 
field lines obtained by adding or subtracting the spin up 
and down field lines drawn in Fig. |SJb). Loop L c winds 
around charge lines but in total it does not wind around 
spin lines and thus it corresponds to pumping of unpo- 
larized charge from left to right. In contrast, loop L s 
winds around spin magnetic lines only and corresponds 
to pumping of pure spin. By pumping spin we mean 
that electrons with opposite spins are transported in op- 
posite directions through the pump. Putting all field 
lines together we reconstruct B° . Assuming it is con- 
tinuous, the conditions for the existence of a finite spin 
current without any charge transport (for infinitesimal 
S) are dS ■ Bf ^ and dS _L Bf. 

Specializing to the case x = a, the effective magnetic 
field in Eq. (|3I|I becomes 

Bf = (0,-d z Q rT ,d y Q (T ). (32) 

The fact that B _L x ensures that there is no pumping 
if a is kept constant. 

We have completed our general classical analysis of the 
spin pumps. In the next section we will apply Eqs. 131|) 
and (|32|) to a particular system of experimental interest. 

V. APPLICATION TO TWO DIMENSIONAL 
ELECTRON GAS 

In the preceding sections we have performed a general 
classical analysis of spin and charge pumps. We have 
argued that with sufficiently strong dephasing our anal- 
ysis is valid. In this section we will apply the general 
theory [Eqs. (O and (J22J] to a 2DEG and propose a 
realization of a spin battery 2 - The pumping parameters 
will be x = a determining the asymmetry between the 
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(a) A 




FIG. 5: (a) Parameter space, (b) Illustrative effective mag- 
netic field lines for spin up and down. The loops Lj and Lj. 
correspond to pumping of spin up and down electrons, respec- 
tively, from left to right, (c) The field lines for charge pumping 
are obtained by adding the spin up and spin down field lines. 
The loop L c corresponds to pumping of charge, i.e., unpolar- 
ized electrons, (d) The effective magnetic lines for spin. The 
loop L s winds around the spin lines but, in total, it does not 
wind around the charge lines. Thus it corresponds to pure 
spin pumping without charge transfer. 



contacts connecting the 2DEG to the left and right leads, 
y = V g e being a plunger gate and z = h = g^sB being 
the Zeeman energy associated with an in-plane magnetic 
field. 

To find B° s we have to calculate the dependence of the 
charge of the 2DEG on the pumping parameters, VQ a 
[see Eq. (|3ip]. The density of particles of each spin in 
the system, n ai is determined by the grand canonical 
ensemble average, and thus it is a function of the chemical 
potentials and of the parameters r — (x, y, z, ...): n a — 
n a (fJ<-\, Mj., r) or zv = /i CT (n|, n^,r). In the 2D case Q a = 
n a eA where A is the 2DEG area. The differentials df, 
dn a for which the system remains at equilibrium satisfy 
d V<j = dl-ia/dn a f\ n _,^ ■ dn a > + V/i CT | nt ,ni ■ dr = 0. 
This leads to the equality 



VQ 



dQp 

dxi 



(33) 



where the thermodynamical density of states tensor 
(DOS) was introduced, 



D„„, = 



(34) 



The quantity (V/v)? is referred as the r$ inverse com- 
pressibility of spin a. 

When the energy is given explicitly as function of n a 
and Ti we can obtain the chemical potential, the r% inverse 
compressibility of spin a and the DOS tensor: 



/'a 

F / 

* rr.rr' 



D = F 



(35) 



Assuming that the pumping parameters influence the 
energy of the system only via the terms — e(rif + n^)V g — 
h(n-\ —rii) then from Eq. Q35JI it follows (V/i CT | nTi „ i ) !/ = 
-1 and (VjUgUt ,nj )z = o~ with a = ±1. Using = 
Dji and Eq. J3lJ we get 



Bf 
Bf 
Bf 



eA(0, a(p a s - D aa ), D arJ + D aW ) 
eA(0,D n -D n ,D n +D u +2D n ) 
eA(0,2D n ~D n 



D 



(36) 



with x — a, y = eV gi z = h and a = —a. Together with 
Eq. I|31l) we can express the pumped current, I a , in terms 
of the DOS tensor. 

There are some interesting relations concerning the 
DOS tensor: in a pumping process in which h (Zeeman 
energy) remains constant we find that only the third com- 
ponent of B cS contributes. Using Eq. (|36|) . the ratio be- 
tween the spin and charge currents in this case is 



B 



off 
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2D 



II 



(37) 



In equilibrium we have /Xf = Hi = /i. MacDon- 
ald showed^ that the inverse magnetic compressibility, 
is given by minus the expression in Eq. (|37|l . 



dh l«T+ ra l ' 



Similarly, at constant V g one can verify that y- — 

^|„ T _ ni . Therefore measurement of charge and spin 
pumping currents reveals thermodynamical properties of 
the system. 

Now we will evaluate the spin and charge pumped cur- 
rents in a 2DEG relying on the Hartree-Fock approxima- 
tion for its energy^ 



E(n h ni,y,z) = 



if+nf e 2 (n T +n x ) 2 



2D 



2C 



8e 2 (n^ 



3/2 



3V^ 



vi n ] + n l) - z ( n ] ~ n l), (38) 



where D n 



2irh 2 /m*. The first term is the kinetic 



energy, the second and third terms are the Hartree and 
the exchange interactions respectively, the forth term de- 
scribes interaction with an external gate potential eV g = 
y and the last term describes interaction with an external 
magnetic field h = g^sB = z. As we will see, the neg- 
ative sign of the exchange term increases the magnetic 
susceptibility at low densities. 



A. Non-interacting electrons 

To discuss the noninteracting case we disregard the 
terms cx e 2 in Eq. (|38|l . One finds then using Eq. Ij35(l 
D 
D 

tive magnetic fields, in units of eADo (charge per unit 



that D 



Equation l|36(l gives the effec- 
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energy), 



B 



off 



(0,-<r,l) 



B] 



oil 



(0,0,2) Bf= (0,-2,0). 

(39) 

The facts that Bf\\z and Bf\\y (y = h, z = V g ) mean 
that a small loop pointing in the z direction (changing a 
and V g ) produces only charge current, and a small loop 
pointing in the y direction (changing a and h) produces 
a pure spin current. (The direction of a infinitesimal loop 
is normal to the plane of the loop.) 

To estimate the pumped currents for oscillating a and 
V g or magnetic field B we use the expression I c<s = 
2T^ 1 eAD 5r c ^Sa derived from Eqs. (JSU and Here 
the energy Sr c is SV g e, when only charge is pumped, and 
Sr s is Sh, when only spin is pumped. 

If we assume that the left-right resistors of the pump 
oscillate with maximal amplitude, 5a = 1, and we con- 
sider a GaAs sample (g=-0.44) of area A = lpm 2 , then 
the charge current obtained for 8V g — 3mV and fre- 
quency t -1 = 10GHz is 1.3/iA At high frequencies it is 
difficult experimentally to produce oscillatory magnetic 
field with a sizable Zeeman energy. Thus, for a magnetic 
field of order SB — lmT and frequency r _1 = lOKHz 
the spin pumped current is very small ~ .5 10~ 15 A. 



B. Capacitive interaction 

The simplest way to consider the Coulomb interaction 
is performed by adding a capacitive term to the Hamilto- 
nian. Thus, adding the second term in Eq. (|38[) one has 
in units of eADo 

Bf = (0,-a,l/(l + 2D Q e 2 /C)) 

Bf = (0,0, 2/(1 + 2D e 2 /C)) 

Bf = (0,-2,0). (40) 

The observation that the magnitude of Bf nas decreased 
reflects the fact that the energy needed to add charge to 
the system is now larger. On the other hand, Bf and 
hence the spin current are unaffected by the capacitance 
term. 

The capacitance between the 2DEG and the gate sep- 
arated one from the other by a distance d is C = 
Then, the reduction factor l + 2Doe 2 jC can be written as 
1 + where af is the effective Bohr radius « 100A in 

GaAs. For a separation of d = 1000 A the reduction 
factor is 20. An estimate for the pumped charge in the 
presence of the capacitance with SV g = 30m V , under the 
conditions specified above, yields I c ~ 0.66/jA 

In contrast to biased transport, where the conductivity 
is proportional to the aspect ratio in two dimensions, 
the current of the pump is proportional to the area A 
as can be seen in the pre- factor of B cS . There are two 
reasons to bound the area of the pump from above. One 
is that spin-flip events may spoil the spin polarization 
in too big samples. This restricts the typical lengths to 



~ jumS The other restriction comes from the adiabatic 
condition, which for the simple circuit of Fig. reads 
RC <C t, where the capacitance is proportional to the 
area A. Rewriting the energy per unit area in terms of 
density, n = + and magnetization, m = n-j- — nj,, as 



E(n, m) 



[to - m (/i)] 
[n - n (V g )} 



1 



4A> 
2{ 1 



V4A) 2(7 



(41) 



where no and too are the average density and magneti- 
zation respectively, we can read off the effective capac- 
itances corresponding to charge, C c = (^ + 2 AD e i 
and to spin, C s = 2ADoe 2 , where C — CA. We see that 
C s > C c , thus the adiabatic condition i?max{C c , C s } <§; 

off 



t becomes A <C 



tions with R e 
A < 400pm 2 . 



eR 

Rl II Rr 



For the above mentioned condi- 
\KVt, the adiabatic bound is 



C. Exchange interaction 



The next step is to include the exchange term, the 
third in Eq. H38|) . Now the DOS tensor depends on elec- 
tron density: F = D and F acr < — S aa <D Q + ^ — 



2e 2 rC 1 / 2 



which is different from the noninteracting 
where the DOS is constant. The depen- 



2D systems, 

dence on density leads to a spontaneous magnetization 
at low densities corresponding to r s = 2,— The spon- 
taneous magnetization occurs in our model due to the 
following facts: (i) at constant density n = rij + nj, the 
exchange term is minimized for maximal \m\, i.e., for full 
polarization; (ii) in sufficiently low densities the exchange 
term dominates, since it behaves as n 3 / 2 and the others 
as n 2 . 

Before the spontaneous magnetization occurs applica- 
tion of magnetic field strongly polarizes the spins (ti| ^ 



7i j) and hence Djf — cx 



1/2 



-1/2 



^ 0. This im- 



balance between the spin up and spin down populations 
at weak magnetic field becomes stronger at low densi- 
ties. Using Eq. I|36|l we see that except for special points 
in parameter space it is not true that B° \\z and Bf\\y. 
From the last fact it follows that one can achieve spin 
current without varying magnetic field. This spin cur- 
rent is proportional to — and becomes larger at 
low densities. 

Following the last observation we calculated I s and I c 
as function of r s = at a constant magnetic field 

of IT. The conditions were similar to those specified in 
Sec. ESI with d = 1000A, SVg = 30ml/ and 6a = 1. 
The results are plotted in Fig. [fjj The charge current 
is almost invariant I c ~ 0.66fiA. On the other hand, as 
explained in the previous paragraph, / s grows rapidly as 
the density is decreased until the ground state becomes 
ferromagnetic at r, ~ 1.4 (the critical r s is now smaller 
due to the magnetic field). 
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FIG. 6: Spin current (big dots) (cc — D±±) and charge cur- 
rent (small dots) (oc 4-fj.j. +2D^i) as function of density 
in a 2DEG according to the Hartree Fock approximation. 



We do not expect the Hartree-Fock approximation to 
be reliable for r s > 1. Monte Carlo calculations^ show 
that, in reality this instability does not occur until much 
lower densities are reached. For r s = 3/4 we have I s — 
15.6nA It should be emphasized that the existence of 
I s is a consequence of the interactions between electrons 
in two dimensions. Notice that this spin current flows in 
the background of a large charge current of 0.64/iA and 
may be hard to detect. 

This difficulty can be resolved by a manipulation of 
the puming contour, using the fact that though the spin 
current is small, i.e., _B| ff <C B^ s , it is rapidly varying 
at low densities, -^Bf > -^Bf. (In Fig. ED one 

can see how these quantities depend on density.) Ac- 
cording to our topological interpretation, a small loop 
in parameter space probes the effective magnetic field 
lines passing through it. In order to probe the deriva- 
tive of B eS , and not its magnitude, we will use a figure 
8 shaped trajectory whose equation in the (a, V g ) plane 
is a(t) = Sas'm(2ujt), V g (t) — SV g s'm(ujt). Basically, this 
trajectory is composed of two adjacent loops with op- 
posite directions. Applying such a pumping trajectory 
for the case of the 2DEG, means the following: in the 
first half of the period (one loop) we pump a charge with 
small polarization, as we found in our previous calcula- 
tion at r s — 3/4. In the second half we pump a nearly 
equal amount of charge to the opposite direction, since 
the direction of the second loop is reversed. In the second 
part, however, the average density is different, thus the 
spin polarization is different. The excess spins may lead 
to a pure spin current. 

The calculation for the pumped currents with the fig- 
ure 8 shaped loop, as defined above with 6a = 1/2 and 
SV g = 50mV yields I c = 1.2nA, I s = 5.2nA, £ = 4.3. 
In order to cancel the residual 1.2nA charge current we 
will consider an asymmetric figure 8 shaped loop. This 
can be done using the parametrization in the (a, V g ) 
plane given above with an additional term in V g equal 
to eSV g \ sm(uit)\, where the asymmetry of the trajectory 
is determined by the parameter |e| < 1. We find that 
5* diverges at Eq = — 10~ 3 . However, getting exactly to 



the point where the charge current is cancelled requires 
a fine tuning of e. For example, to obtain > 40, the 
parameter e should be close to Eq as \e — eq\ < 2 x 10~ 4 . 
This means that a resolution of about 10[j,V in the gates 
voltage is required. 



D. Temperature effects 

We would like to estimate the temperature needed for 
the operation of the spin pump discussed in Sec. IV CI 
Consider the quadratic part of Eq. I|38|) . given in Eq. 141|) . 
The thermal smearing of the magnetization m is equal to 

ArriT — \j DaK A^ - (The presence of the exchange term 

does not change Amy significantly) . In the pumping pro- 
cess the parameters oscillate and change the magnetiza- 
tion m during the period with an amplitude Amp, which 
is related to the pumped spin current by Amp ~ 
We impose that the thermal smearing of m is small rel- 
ative to the parameters induced oscillations in m, i.e., 
Am T < Amp. This yields T < i 1 ^) 2 Kb d oA = 
where I s — 5.2nA was used. Finally, the thermal smear- 
ing of the density n is smaller than Amy by the factor 

= 4.4 for the above specified conditions. 



VI. SUMMARY 

We analyzed classical non-coherent pumps and dis- 
cussed their connection with quantum pumps using a 
dephasing model. Our classical method can include 
electron-electron interactions and allows to study the ef- 
fect of exchange interaction on spin pumping. We ex- 
pressed the classical pumped spin current in terms of the 
thermodynamic DOS tensor of the system, and gave a 
topological interpretation to it in terms of effective "mag- 
netic" flux through trajectory loops in the space of the 
parameters that control the pump [see Eq. ©]. We an- 
alyzed in detail the case of 2DEG GaAs and found that 
any combination of charge and spin currents can be ob- 
tained by choosing appropriate trajectory in parameter 
space. In particular one can choose a trajectory that cor- 
responds to a pure spin current, which has magnitude of 
order of nano-Ampers. 
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APPENDIX A: PUMPED CURRENTS VS. BIAS 
CURRENTS 

According to Einstein's relation the current of spin 
o =T, I is 

4ias, CT /Vbias = e D ^~g^ x I as P ect ratio], (Al) 

where -Djf is the diffusion constant. Using the definition 
of the DOS tensor, Eq. I|34|l . we find that in equilibrium 
(Mt = Ml = M) we have = D aa , . 

Assuming that the diffusion constant is spin indepen- 
dent we find, [similar to Eq. (|37|l for the case of pumps] 
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(A2) 



This shows that bias currents and classical pumped cur- 
rents have similar expressions in terms of the DOS tensor, 
cf Eq. ■ Now we will show that a pure spin current, on 
which we elaborated in Sec. IV Cl for the case of pumping, 
can be obtained by an oscillatory bias together with an 
in-phase oscillating gate. 

Consider an AC biased 2DEG in a constant magnetic 
field, where a gate V g oscillates in phase with the bias: 
VbiasW = V b sva(u)t), VgateC*) = V g sm(ujt). Due to the 
oscillating gate a rectified DC component in the current 
will appear. Using Einstein's relation, this DC current is 
given by 

(I c ) = e 2 D i{ ((D n + D u +2D u )(t)V bias (t)) x [a. r.] 
(I s ) = e 2 D if {(D n -D u )(t)V hia *(t)) x [a- r.] (A3) 

where () denotes time average over a period and we as- 
sumed that the diffusion constant does not depend on 
Vbias or V g 

ate • 

Near the magnetic instability the spin current oc — 
Dn is very susceptible to changes in V a while the charge 
current cx D^+Dn+2D^ is nearly constant (see Fig.|HJ). 
As a result, the charge flowing in one half of the period 
cancels the charge flowing in the second half, while a 
finite amount of spins is transferred from left to tight. 

This rectification 1 — effect is based on the same mech- 
anism as the adiabatic pumping effect considered in 
Sec. IV CI Both effects are proportional to the interac- 
tion strength. 

We will show now that the two effects of pumping and 
rectification (bias) produce currents of the same order, 
provided that (i) the pumping frequency saturates the 
adiabatic condition, r = RC and (ii) the pumping gate, 
the rectification gate and the bias oscillate with similar 
amplitudes. 



First consider a simple pumping loop in which the 
asymmetry parameter a and the pumping gate oscil- 
late out of phase with 6a = 1, 6V g = V. When the 
adiabatic condition is saturated, r = RC , the pumped 
current in the circuit of Fig. given in Eq. i|12|l. is 
J pump = T~ 1 CVSa = V/R. Using Ohm's law we see that 
the DC current flowing due to a DC bias V in the same 
circuit is of the same order, i.e., /bias — V/2R, where we 
assumed R = Rl || Rr — Rl/2 — Rr/2 for simplicity. 

Second, let us compare the spin currents for pumping 
with a figure 8 shaped loop and for rectification (i.e., 
bias voltage) with an oscillating gate. We can apply the 
above argument for each half period of either pumping 
or rectification processes: for pumping, each half period 
consists of one of the two loops of the figure 8. In each 
such loop the amplitude of the gate oscillation is half 
of the oscillation in the entire period. For the case of 
rectification, in each half period the averaged bias is of 
the order of the maximal AC amplitude. In both cases, 
the density varies between the two halves of the periods: 
the voltage providing different densities corresponds to 
(i) the distance between the centers of the two loops in 
figure 8 for the case of pumping; (ii) the oscillations of the 
gate in the case of rectification. Thus when the pumping 
gate, the rectification gate and the bias oscillate with 
similar amplitudes similar spin currents will be created 
in the pumping and in the rectification processes. 

To conclude, unless biasing of the system is unwanted, 
one can obtain spin current either by the pumping ef- 
fect discussed in Sec. IV CI or by the rectification effect 
combined with an oscillating gate described here. 

On the other hand, as is seen in Appendix [BJ for long 
systems adiabatic pumping is more feasible than appli- 
cation of bias, since it requires only small voltages. 



APPENDIX B: CONNECTION WITH 
NON-LINEAR RESPONSE THEORY 

In this section we compare our results to other 
theories22i2i that evaluate non linear response to an ex- 
ternal potential at wave number k and frequency w. 

We start by analyzing the circuit of Fig. [3] in the con- 
tinuous limit, with the assumptions (i) R m — R(V g m ) 
and C m = C(V g rn ) so that the varying parameters are 
{V g m }; (ii) periodic boundary conditions: R\ is con- 
nected to Rn- 

The continuous version of the current conservation in 
the junctions of Fig.^is, in the adiabatic limit, 



dl(r, t) 
dr 



= -C(V g (r,t))V g (r,t), u;RCk- 2 ^l, (Bl) 



where I(r) is the current along the circuit and C, R are 
the capacitance and resistance per unit length respec- 
tively. To estimate the adiabatic condition we use the 
characteristic length scale The periodic boundary 
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conditions imply 



Eq. we get the charge per period for small area §: 



drI(r,t)R(V g (r,t)) = 0, 



(B2) 



where I is the length of the system. Integrating Eq. I|B1|) 
and substituting into Eq. i|B2|) we have: 



CR' 
kR Js 



dxdy 



R'Cu 2 



(B7) 



I(0,t)lR(V g (t)) = / drR(V g {r,t)) 



x / dr'C(V g (r',t))V g (r',t),(B3) 
Jo 

where V g (t) = r 1 J ' cW s (r,t). 

We proceed by assuming that the potential V g (r, t) de- 
pends on two parameters x, y as 

V g (r, t) = x(t) sin(fcr) + y(t) cos(fcr), (B4) 

and take C(V g (r, t)) — C to be a constant. The integral 
over r' in Eq. I|B3|I yields 

C 

— [x(t) (1 - cos(£;r)) + y(t) sin(Ar)] , (B5) 
and the integral over r yields 

HO, t) = H (x(*)l/(t) - l/(t)i(*)) + (B6) 

where we expanded i? around = as R(V g ) = R + 
R'Vg. Integrating Ea. ljBBj) over one period as done in 



in the last equality we have assumed that x(t) — 
V2Vcos(ujt), y(t) = V2Vam{ujt). 

We can compare our results to the non liner re- 
sponse theory. Ignoring coherent effects^ the current 
response to a time and space dependent potential V(r, t) 
is j(r,t) = —cr{n(r,t)}\7V(r,t). If V(r,t) induces a den- 
sity polarization given by Sn(k,uj) — — Tl(k, u>)eV(k, to), 
the DC response current is given by22*2i 

k,u 



We take now the diffusive polarization operator— in 
the limit D i{ k 2 > w, this gives Sffl(fc,w) - - " 
Thus, 



D it k' 2 dn ' 



__ _da_dfL_l_ 2dm, u_ v 2 

JDC dn den D i{ c d^ k v 

_ 

R 



c ±v\ 



(B9) 



e z dn/ /. 



which is identical to Eq. I|B7(I . Since D-^ = 

\/{RC) the limit -Diffc 2 3> u> is identical to the adiabatic 
limit, Eq. (JbTJ. 
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